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Abstract 

The concept of a propagator is useful and is a well-known object in diffusion NMR 
experiments. Here, we investigate the related concept; the propagator for the magnetisation 
or the Green’s function of the Torrey-Bloch equations. The magnetisation propagator is 
constructed by defining functions such as the Hamiltonian and Lagrangian and using these to 
define a path integral. It is shown that the equations of motion derived from the Lagrangian 
produce complex-valued trajectories (classical paths) and it is conjectured that the end¬ 
points of these trajectories are real-valued. The complex nature of the trajectories also 
suggests that the spin degrees of freedom are also encoded into the trajectories and this idea 
is explored by explicitly modeling the spin or precessing magnetisation by anticommuting 
Grassmann variables. A pseudoclassical Lagrangian is constructed by combining the diffusive 
(bosonic) Lagrangian with the Grassmann (fermionic) Lagrangian, and performing the path 
integral over the Grassmann variables recovers the original Lagrangian that was used in the 
construction of the propagator for the magnetisation. The trajectories of the pseudoclassical 
model also provide some insight into the nature of the end-points. 
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1 Introduction 

Nuclear magnetic resonance (NMR) is a versatile experimental technique that has acquired a 
vast and varied collection of methods to measure what might broadly be categorised as molecular 
composition, structure, and dynamics. Within the last category are the translational motions 
of diffusing and flowing particles, and it is with these types of motion that we will be concerned 
with here. In particular, we will focus on the NMR measurement of translational motions of 
fluid particles mm by methods employing field gradients. 

A common method of modeling NMR experiments of translational motion is via the prop¬ 
agator [3]. The propagator is otherwise known as the conditional probability density [1] and is 
obtained as the Green’s function to the underlying differential equation of the motion, usually a 
diffusion or Fokker-Planck equation [S] of some form. Once the propagator is determined it can 
be used to average the NMR phase factor exp(i f dtjB(x,t)) to obtain the time-dependence of 
the macroscopic magnetisation due to motion in the magnetic field 

In the context of diffusion or the Fokker-Planck equation, the propagator is an object which 
describes statistics of the stochastic process. For Markovian processes [1], which will be assumed 
henceforth, the propagator together with the initial probability (marginal) density contain all 
the statistical information available. The propagator is thus an important function and there are 
many approaches available for its calculation. One such method is via the powerful technique 
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of path integration [SHH] which was originally created to calculate the analogous quantity in 
quantum mechanics [MU- In both these contexts, path integration is widely used, as it is in 
related subjects such as statistical physics. It is therefore surprising that it has been very rarely 
exploited in the context of NMR, with the exception of Le Doussal and Sen m- 

There are possibly a few reasons for this absence of path integral methods. The first reason 
is probably because the propagator of the diffusion or Fokker-Planck equation is often more 
easily derivable by other methods in cases where the propagator is non-Gaussian, such as in 
the case of restricted or reflected diffusion m. where an eigenfunction expansion, for example, 
works well. Path integral techniques deliver exact solutions in the case of Gaussian propagators 
but most of those cases are already known and derivable by other methods—for example, by 
cumulant expansions m- Therefore, in regard to the propagator for the underlying stochastic 
process, unless there are any benefits of the perturbation methods particular to path integration, 
there might not be any circumstances in which path integration is advantageous in the context 
of NMR. 

Suppose, however, that instead of averaging the NMR phase factor with the stochastic 
propagator, the propagator for the magnetisation is calculated directly. This is essentially what 
was done by Le Doussal and Sen [12] who calculated the magnetisation propagator for diffusion 
in a time-independent parabolic magnetic field by analytically continuing the analogous quantum 
mechanical propagator for a particle in a quadratic potential mm- 

It may seem that there is no significant difference between these two approaches, besides the 
fact that the propagator in the latter case becomes a complex function, but there is a subtlety 
which is not apparent until a path integral is formulated directly from the partial differential 
equation of the magnetisation. In the simplest case of diffusion in a linear field gradient, the 
term iyg • xM(x,t) is added to the diffusion equation for the magnetisation M(x,t) to form 
the Torrey-Bloch equation. It is this term which couples the diffusive motion with the spin 
(magnetic) degrees of freedom and introduces a non-trivial element to the dynamics. 

This non-trivial element is not apparent until one examines the analytical mechanics and the 
ensuing equations of motion via the Hamiltonian or Lagrangian. Usually, whether one considers 
classical or quantum mechanics, motion in a linear potential (in the absence of any damping 
terms such as viscosity) leads to accelerated motion. This is also true for diffusion NMR but 
it is found that the acceleration is imaginary and thus the “classical” or “particle” trajectory 
given by the Hamiltonian or Lagrangian dynamics is generally complex-valued. 

This, at first, seems like an inconsistency: How can the position and velocity variables be 
complex? The answer is surely that we have unwittingly also encoded the spin degrees of freedom 
into the variables and this, in turn, leads to the intriguing possibility that we might also be able 
to include the spin degrees of freedom explicitly in the propagator. The results presented here 
are a preliminary investigation of this idea and related issues. 

In section[2|we begin by demonstrating the construction of a propagator for the magnetisation 
in a simple case of a diffusion NMR experiment. The propagator is constructed as a path integral 
over (paths of) the ordinary position variables but it is shown that the classical trajectory—the 
extremum path of the action—is complex-valued, and it is then conjectured that the end-points 
of these trajectories should be real-valued. 

Section [3| introduces the anticommuting Grassmann variables and shows that a simple 
Lagrangian produces equations of motion which mimic processing magnetisation. A pseudo- 
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classical Lagrangian (section 0]) is then constructed from both the ordinary position variables 
and the Grassmann variables. An analysis of the equations of motion from this Lagrangian 
then lend some support to the notion that the end-points of the complex trajectories should be 
real. In sections 14.1H4.3I it is shown that by performing the path integration over the Grassmann 
variables the original propagator (involving just the position variables) is recovered. 

In the final discussion, we mention some ways of extending the results. In particular, we 
discuss the possibility of incorporating radiofrequency fields into the propagator and—one of 
the main motivating factors behind this work—the possibility of constructing a supersymmetric 
model. 


2 Path integral for an NMR propagator 


2.1 Torrey-Bloch equations 


The Bloch equations [ni[E] describe the macroscopic magnetisation M(t) of an ensemble of 
spins in an applied magnetic field B(t) 


dM^ 

dt 

dMy 

dt 

dM^ 

dt 


7(M X B)^ 
7(M X B)j^ 
7(MxB), 


T2 ’ 

My 
T2 ’ 

M, -Mq 
Ti 


( 2 . 1 ) 


where 7 is the gyromagnetic ratio, Ti and T 2 are the spin-lattice and spin-spin relaxation times 
respectively, and Mq is the equilibrium magnetization. In the absence of radiofrequency pulses, 
B(t) = B{t)z is purely polarized in the z-direction but might vary temporally (and spatially). 
In this case, the z-component of the magnetization is uninteresting (as far as our particular 
interest goes) and we can focus solely on the x and y-components, so that 


^ (MA 

dt \My) 





^ (MA 
T2 \MyJ 


( 2 . 2 ) 


where 



(2.3) 


which is the generator of the Lie group SO(3), reduced to the x,y components, i.e. the 
generator of rotations about the z-axis. 

More usually, a complex magnetisation is defined such that M± = ± iMy and we can 

then transform the above vector equation into two decoupled components such that 


where 


dt \M-) 





1 [M. 


T2 \M. 


1+ 



(2.4) 
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is the su(2) Pauli matrix corresponding to the so(3) matrix g^. 

To generalize the differential equation (j2.4p to include diffusive and coherent motion we can 
combine the above with the Fokker-Planck equation [5] 

f)P 

— = V2(DP)-V(vP), (2.6) 

where P is a density distribution and, in the most general case, D is a temporally and spatially- 
dependent diffusion tensor, and v is a temporally and spatially-dependent velocity field. 

We shall consider the less general case in which D = P is a temporally and spatially- 
independent scalar, in which case the Fokker-Planck equation reduces to the simpler form 

dP 
'dt 


= - V (vP). 


(2.7) 


The right-hand side of this partial differential equation can be combined with, say, the M+ 
component (which will henceforth be denoted more simply as M) of equation (j2.4p to form the 
Torrey-Bloch equation |141I19] with an advection term 


dM 


= DV'^M - V (vM) - i'yB{t)M - 


M 


( 2 . 8 ) 


dt ' Ts’ 

where it should also be mentioned that T 2 could, in principle, be spatially and temporally 
dependent. 

Let us make one further simplification. The applied magnetic field B{t) is composed of two 
components. The first is a large, time-independent, spatially-uniform magnetic field Bq. This 
is the field generated by, for example, a superconducting magnet which polarizes the individual 
spin magnetic moments and leads to a macroscopic magnetisation within the sample. During a 
typical diffusion measurement, linear magnetic field gradients are applied at various times which 
encode the motion of the sample molecules. These gradients may vary the ^;-component of the 
magnetic field linearly in any directioE0 x with strength g. Thus the total magnetic field is given 
by 

P(x,t) = Po + g(i)-x. (2.9) 

The constant Bq can be neglected, however. This term merely modulates M with the oscillating 
factor exp(—iyPo^) and is not important in the dynamics that we will investigate below. Ignoring 
the Bq, we arrive at our defining equation for the magnetization M, 


dM 

~dt 


DV'^M - V - i-fgit) • xM - 


( 2 . 10 ) 


In the next two sections, for simplicity, we will focus on the diffusive part of (j2.inh . setting 
V = 0 and R 2 = l/r 2 = 0. Thus, 


dM 

~df 


= DS/^M — i^g{t) ■ xM 


( 2 . 11 ) 


will be our starting point. However, a more general differential equation is considered in Ap¬ 
pendix |Al 


^Here x denotes a general position vector. In the following sections, index notation will be used for convenience. 
There, the components of x will be denoted by a:“, a = 1, 2, 3. This is equivalent to x = {x, y, z)'^. 
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2.2 Hamiltonian, Lagrangian, equations of motion 

Clearly, the PDE for M, equation (I2.1ip is not really an equation of energy. Nonetheless, it is 
analogous (to some extent) to a Schrodinger equation and quantities such as the Hamiltonian 
and momentum operators may be defined analogously. 

To begin, we define the momentum operator a^ 

so that our phase-space coordinates are with a = 1,2,3. The Hamiltonian, correspond¬ 

ing to the operator —d/dt, is then 


H{x,p) = Dpap^- + igaX° 


(2.13) 


where index notation (and Einstein summation convention) is now being used, and where 7 has 
been absorbed into the definition of g. 

To define the Lagrangian we need an expression for the velocity This is found using one 
of Hamilton’s equations of motion (see Appendix]^, i.e. 


dH 


Xa = 


dH 

= *77-7 = *2T)pa- 


d{—ip°‘) dp 

The Lagrangian is then defined via the Legendre transform 


(2.14) 


L{x,x) = —ipaX°' — H{x,p) 
= -^XaX^- -igaX° 


(2.15) 


x) 


Lagrange’s equations of motion satisfy (in the absence of higher order time-derivatives than 


dL d dL 


dx°- dt dx°‘ 


= 0 , 


(2.16) 


from which we obtain 


= 




i:“ = + i2D 


rdV(f'), 

Jto 


= 


Xq + v°‘{t — to) + i2D f dt' f dt"g°'{t"), 


(2.17) 

(2.18) 

(2.19) 


'to 


'to 


^As mentioned in the text, the diffusion or Fokker-Planck equation is a PDE for a stochastic process, not 
a quantum process, and therefore the usual definitions of quantities such as the Hamiltonian, Lagrangian, and 
momentum are slightly different. See Appendix]^ for a more general derivation of the Lagrangian, Hamiltonian 
and associated definitions. 
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where Xq and f “ are constants of integration. In the case of constant g, the equations of motion 
simplify to (setting to = 0) 


i2Dg\ 

(2.20) 

-h i2Dg°‘t, 

(2.21) 

Xq + v°‘t + iDg°'t^, 

(2.22) 


where it is clearer that equation (|2.22h for the trajectory in a constant magnetic held gradient 
has an imaginary acceleration and is, in general, complex-valued. Note that Xq and are 
also generally complex numbers. By separating the real and imaginary parts for one component 
(one dimension), such that x = Re(x“), y = Im(x“) for hxed a, the trajectory can be made 
analogous to the two-dimensional parabolic motion of a body in a linear potential 


X = xq + vt, (2.23) 

y = yQ + ut + Dgt^. (2.24) 


The form of y is also the same as that of the NMR phase as a particle moves at constant 
velocity v in the direction of the gradient, 

4> = </>o + / dt'gx{t'), 

Jo 

= (l >0 + gxot + ^gvt^, (2.25) 

but there appears to be no obvious link between the imaginary component y and the phase 4> 
and, therefore, the imaginary part of the trajectory is not simply proportional to the phase. 


2.3 Path integral, propagator, magnetisation 

With the Lagrangian defined in (|2.15l) . the propagator (or Green’s function) for the magnetisa¬ 
tion of equation (12.lip can be calculated using the technique of path integration [THTO] . Denoting 
the propagator by G{xt,T; xo,to), 


G{xt,T;xo, 



Vx , 


(2.26) 


where f T>x represents an inhnite-dimensional integral over all paths x{t) beginning at xo,to and 
ending at xt, T. The action S is given by the integral of the Lagrangian; 

S= f dtL{x,x). (2.27) 

Jto 

Since the Lagrangian (I2.15p is quadratic in x and x, the path integral can be calculated 
exactly. The action is expanded in a variational manner 


S = So + 6S + 6^S, 


(2.28) 
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where, because the Lagrangian is quadratic the variation terminates at the <5^ term. The hrst 
term Sq is the action evaluated along the path given by the equations of motion, equation (12.191) . 
In the context of quantum mechanics this would be the classical trajectory. By design, the 
first variation 6S = 0 due to the Euler-Lagrange equation (12.161) . Finally, the third term 
6‘^S is independent of x and x and corresponds to the contribution to the path integral of 
the fluctuations around the classical trajectory. This term produces the normalisation to the 
propagator, 

G {xT, T ; xo, to) = i^' dxTdx^ T, xq, to)), (2.29) 

and is the square-root of the Van Vleck determinant m of the matrix d'^So/dx^dxQ. 

Using the trajectory, equations (|2.18l) and (I2.19j) . Sq can be evaluated as 


5o 


{xT - Xpf 
AD{T - to) 

ixT [ dtgait) + i^^—[ dt f dt'ga{t') 

Jto -L — to JtQ Jto 

D Tdt I fdt'g{t')\ Tdt fdt'g{A 

Jto KJto ) — to \^Jto Jto ) 


and the pre-exponential normalisation factor in (|2.29l) is 


(2.30) 



1 9^50 \ 

27r dxxdxQ J 


1 

{47rD{T-to)f^' 


(2.31) 


Note that in most NMR diffusion experiments the effective time-dependence of the gradient is 
chosen such that its integral at the time of measurement is zero. Thus, usually the second term 
on the right-hand side can be neglected. 

In the case that g is constant, the propagator is given by 


G {xT,T;xo,to) = 


■ ^4D(T-°tl) -^('^-^o)9a{xT+Xo)°--^g^{T-tof 


(AirDiT - to))"/" 


(2.32) 


which compares analogously with the quantum mechanical propagator for a particle in a constant 
linear electric field [am] and more exactly with the results of Stoller et al. m and Le Doussal 
and Sen [T2|. However, the position variables in these references are real whereas, apparently, 
our variables x^ and Xq are complex. Thus, we see that if we arbitrarily set the end-points of the 
trajectories to be real numbers (whilst still allowing the rest of the trajectories to be generally 
complex-valued) our propagator agrees with previous results. 

At this stage, it is not clear whether this real constraint on the end-points is necessary and, 
if so, whether it may be specified less arbitrarily than by fiat. Conversely, if it is not necessary, 
it suggests that the imaginary part could be measurable experimentally. For now, however, we 
shall form the following conjecture: 






















Conjecture 2.1 (Reality Constraint). Let the path x(t) : [0,T] be parameterised by t €z 

[0,T]. Then, to obtain a valid propagator for the diffusion NMR experiment, it is a requirement 
that lm{x(0)} = Im{x(T)} = 0. 

To complete the picture before we move on into less well-known territory, let us compute 
the attenuation of magnetisation that occurs in this case; an NMR diffusion experiment for free 
diffusion in a linear, time-dependent, field gradient (which includes a time-dependence of the 
direction). For convenience, and as is usually the case, we will assume that the time-integral of 
the gradient over the time T from initial excitation, to = 0, to measurement, t = T, is zero. The 
magnetisation at t = T is then given by 

/ oo roo 

d^XT / d^xoG {xt,T;xo, 0) M{xo,0). (2.33) 

-OO J —oo 

With the time-integral constraint on g(t), the action given above in equation (12.3011 is only 
dependent on the difference xt — xq. Thus, by changing variables from xt to r = xt — xq, the 
two integrals can be separated; 

/ OO POO 

dVG(r,T) / d\oM(a:o,0). (2.34) 

-OO J —OO 


Letting the integral over xq be denoted by Mq (which may be complex), the integrations over r 
can be readily performed (being merely Gaussian integrals) to give 


(M(r)) 




(2.35) 


as was derived by Karlicek and Lowe m and initially (in slightly less general form) by Stejskal 
and Tanner [TS] by very different methods (also see Kuchel et al. |19|L 

Thus, we see that the usual expression for the attenuation of the magnetisation in a diffusion 
experiment is derivable from a path integral method, and this provides us with some evidence 
that the method is correct. In fact, propagators can also be derived in more general situations, 
such as with diffusion tensors, velocity fields, potentials, and non-constant relaxation terms (see 
Appendix|A]for a more general method of forming the path integral) but, except in special cases, 
the path integral is often not exactly calculable when the Lagrangian involves higher powers than 
the second. 

Here, however, we shall not be concerned with these more general cases. Instead, the inten¬ 
tion is to investigate the possibility of explicitly including the spin degrees of freedom into the 
path integral in such a way that they account for the imaginary part of the classical trajectory. 


3 Grassmann variables and spin 

The imaginary parts of the trajectories, or classical paths, x°‘(f) given in section [2.21 seem to 
be related somehow to the spin degrees of freedom, or the precessing magnetisation. If this is 
true, we should be able to model the spin components explicitly in some form and, in doing 
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so, this might reveal the nature of the imaginary part of the paths. Non-relativistic spins 
are usually described by Pauli matrices and their su(2) algebra. However, incorporating Pauli 
matrices into path integrals is not straightforward and, besides, it is not clear that modeling the 
spin components quantum-mechanically is the correct way to proceed; after all, the observed 
magnetisation is of a macroscopic nature. 

An alternative is to use Grassmann variables These are often easier to handle within 

path integrals despite their anticommuting properties and can be used to model fermionic fields 
and spin at the classical level |221424] . Scholtz |25] has also shown that it is possible to relate 
path integrals of Grassmann variables to bosonic path integrals over representation variables of 
su(2). Here, we shall consider only Grassmann variables. 

3.1 The Grassmann algebra 

The real Grassmann algebra r(n) of order n is the algebra over M generated from the unit 
element 1 and the n Grassmann variables 0*, i = 1... n, which satisfy the anticommutation 
relations 

{e\e^ = 9^e^ + 9^9^ = 0 . ( 3 . 1 ) 

In particular, the anticommutation means that the variables are nilpotent, 

0*0* = -0*0* = 0, (3.2) 

and this implies that only a finite number, 2**, of monomials are possible; this is the dimension 
of the algebra. 

Specialising to the case n = 3, which is the case we will be interested in, dimr(3) = 8. The 
full set of monomials are 


{!}, {e\9^, 9^}, {0^02, 0203, 0301}, { 010203 } ( 3 . 3 ) 

and they partition the algebra into two parts; even and odd, or (in terms of a mod 2 grading) 
degree 0 and 1 respectively. Thus, r(n) = r(n)Q 0 r(n)j. Specifically, deg 1 = deg 0*0-1 = 0 and 
deg 0 * = deg 0 * 0 - 10 ^ = 1 . 

Further properties of Grassmann variables—particularly in regard to integration—can be 
found in Appendix |Bj 


3.2 A Grassmann Lagrangian 

Now we shall posit a Lagrangian using the Grassmann algebra r(3). Let [81 1^126] 

L(0,0) = 0i0'-H*(^)e^,■fc0l0^ (3.4) 

where H*(t) is a time-dependent vector and €ijk is the antisymmetric Levi-Givita symbol; 


^ijk 


+ 1 , if {i,j, k} is an even permutation of { 1 , 2 ,3} 
< —1, if {z,j, A;} is an odd permutation of (1,2, 3} 
0 , if any indices are repeated. 


(3.5) 
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Using the Euler-Lagrange equation (j2.16p . and being careful with the commutation of odd- 
degree quantities, such as 0 * and djdO'^, we find 

= -eiPB\t)0^. (3.6) 

By expanding the components, it is easy to see that this is informally equivalent to 0 = 6 x B(t) 
and thus the 6 variables precess around the vector B(t) just as with the Bloch equations given 
in 

If we restrict B(t) to be fixed in direction such that B(t) = (0,0, B{t)) the solution to the 
above differential equation is 


9^ = el (3.7) 

'ipit) = -00 exp dTB(r)^ , (3.8) 

where we have introduced the complex odd Grassmann variables ip = 0^ + id‘^ and ip = 9^ — i9'^. 
Following Kleinert [ 8 ], an even-degree quantity can be defined by 

h = -^^ijk9^9^, (3.9) 

and, with this, the interaction part of the Lagrangian above becomes 

B^eijk9^9^ ^ liB^ = 1-^, (3.10) 

which is analogous to the Zeeman interaction of the total spin angular momentum I with a 
magnetic held B. The equations of motion for I corresponding to those of (|3.7p and (|3.8I1 are 

h = = ]^9lel = ^ipoipo, (3.11) 

I+{t) = \9liP{t). (3.12) 

where we have dehned /+ = /i -|- il 2 - Note that whilst the hrst of the above two equations 
looks incorrect, it can be easily verihed by the fact that d( 0 ^ 0 ^)/dt = 0 using the equations of 
motion (13.61) with = B^6\. Thus, is constant and proportional to the amplitudes of the 
precessing motion of 9^ and 0 ^. 

Note that if the Grassmann variables 9^, 9“^, 9^ or the quantities Ii,l 2 ,I^ are regarded as 
vectors or states, they can be rotated amongst themselves by dehning rotation operators 

Ra = (3.13) 

whose commutation relations form an algebra which is homomorphic to the su(2) Lie algebra, 

[Ra, Rb] = Rab Rc- (3.14) 
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4 Pseudoclassical mechanics 

Having established that we can model the spin degrees of freedom (or something akin to the pre- 
cessing magnetisation) via a triplet of (fermionic) Grassmann variables 0*, we now ask whether it 
is possible to combine the Grassmann variables with the (bosonic) spatial variables to obtain 
equations of motion similar to those obtained in the original analysis, equations (I2.17h - (l2.19li . 
or ()2.20l) - (j2.22l) in the case of a time-independent magnetic held gradient. 

We shall attempt to simply combine the Lagrangians for a free particle, L ~ —XaX°‘, with 
the Lagrangian in (13.41) for the Grassmann variables by adding them with an extra parameter 
ia; 

L{x, X, 9,9) = -|- ia9i9^ — iaB'^{x, t)eijk9^9^, (4.1) 

where we have allowed B to be explicitly spatially-dependent so that the two sets of variables 
and 0* will be coupled. 

The equations of motion are 


gfji 

Xa = ia2D-^eijk9^9^, (4.2) 

9m = -eijmB^9K (4.3) 


As previously, let us specify the magnetic held to be polarized in the z-direction, i.e., B = 
(0,0, B^). This hxes the index z = 3 in (14.2D and (14.3D . We then specify B^{x,t) = B{x,t) = 
gi,{t)x^ to be linearly varying in space. Then 


dB^ 

dx°' 



9a6l 


(4.4) 


Inserting this into the equations of motion, we see (as before) that 9^ decouples and we have 


Xa = ia4:Dga9l9l, 
9i' — 9aX €iiji9^ , 

O 3 = 0 , 


(4.5) 

(4.6) 

(4.7) 


where Cj/j/ is the Levi-Givita symbol on two indices such that ei 2 = 1 = —£ 21 ) £11 = £22 = 0, and 
the indices i' and j' are reduced indices taking values 1 or 2 only. We have also used the result 
that 9^9^ = 9q9q, a time-independent quantity. 

The solution for the Grassmann variables is very similar to that in the previous section: 

9^ = 9l (4.8) 

ip{t) = -ipo exp (^-i j drc/a(T)x“(r)^ , (4.9) 

where the exact time-dependence of is dependent on the trajectory for a;“(t), which is given 
by the solution to 

= ia8Dg^{t)l3. (4.10) 
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Comparing this equation of motion for with that of (I2.17jl . it is tempting to suggest that 
a /3 = 1/4. However, this notion overlooks the fact that Is is an even Grassmann parameter. 
Therefore, this suggestion does not make sense unless a can be regarded as an operator of some 
appropriate form. Despite this, if we do entertain the possibility that a and Is can be regarded 
as real numbers, and that the /j represent components of a spin angular momentum vector, then 
for spin -1 we expect that /a should give the component of the spin in the z-direction, and thus 
/3 = 1/2. If this is correct, a = 1/2, which is the usual value of the coupling constant found 
in pseudoclassical [20] and supersymmetric models HlEol. Nonetheless, this subtlety requires 
clarification. Although an answer to this possibly resides in the work of Scholtz [2S|, it will not 
be pursued here. 

Instead, if we take the above equation of motion at face value, it implies that the trajectory for 
develops a component which is not only imaginary but also Grassmann-valued. Following 
Junker and Matthiesen [2711^ . and Junker [2^, we can thus split the trajectory into a so-called 
quasi-classical path x“(f) and a Grassmann-valued part, 

x“(f) = -/fg“(f)/ 3 , (4.11) 


where both and are real functions of time. The quasi-classical path is just the 

trajectory due to the diffusive part of the motion in the absence of any coupling to the spin 
degrees of freedom and is given by the real solutions to = 0 , and corresponds to the trajectory 
in (|2.23l) . 

Recalling the definition of 13 from (jJ.lip . and inserting (14.lip into the trajectory (14.9p for 


= •00 exp -i 


i / dTgair) 
Jo 




= 00 (^1 + ^ J dr5ra(r)g“(r)0o0o^ exp j drffa(r)x“(r) 


= 00 exp -i 


i J dr5ra(r)F(r)^ , 


(4.12) 


where we have used the properties of Grassmann variables to expand the exponential and finally 
eliminate the ( 7 “(f)-term completely using (0o)^ = 0. Thus, we see that the trajectory for 
0 (f) is only dependent on the quasi-classical trajectory x“(f) which, being a real function of 
time, implies that 0 evolves only via its phase, which corresponds exactly to the phase as given 

by urmi . 

The observation that the trajectory for 0(f) and the Grassmann part of the action is inde¬ 
pendent of <?“(f) has a further consequence because we can clearly add an arbitrary function 
A“(f) to Q'“(f) without changing anything (with regard to the Grassmann variables). However, 
there is a constraint on A“(f). From (|4.11h and (|4.10h we also have 


from which A“ = 0, and therefore 


a8Dg^{t)h, 

(4.13) 

■ Ag J- 

(4.14) 

13 




















for arbitrary constants Aq and or, alternatively, constant end-points of the trajectory for 
A“(t): Aq and A^, so that = (A^ — Aq)/T. 

These gauge-like degrees of freedom coincide exactly with the troublesome imaginary parts 
of the trajectory for that was the subject of Conjecture 12.11 and correspond to shifting 

the spatial coordinates q —>■ q -|- iA^Q. However, the propagator (for example, the simple 
propagator of ()2..'12M is generally not invariant under this shift, and neither does it appear to 
change as a result of some physical gauge transformation. Therefore, it appears that in order to 
preserve observables, we must set A^ q = 0 . 

4.1 Grassmann path integration 

The fact that we have almost recovered the original equations of motion (I2.17P for x“ suggests 
that our modeling of the spin degrees of freedom has a chance of being correct or, at least, 
sufficient. However, our work is not yet complete. The equations of motion (I4.10p in the 
previous section are Grassmann-valued and therefore what remains is to perform a path integral 
over the Grassmann variables. 

The path integral over the Grassmann variables is not quite as straightforward muMj as 
with the spatial (commuting) variables. In particular, note that if the equations of motion (|4.6I) 
and (14.71) are substituted in to the Lagrangian (14.11) the Grassmann part of L vanishes. This 
would imply that the propagator of the Grassmann variables is a constant, independent of the 
initial and final values, which cannot be correct. The source of the issue is that the Lagrangian 
is somewhat ambiguous and despite appearing to be formulated in terms of configuration-space 
variables 0“, it is really a phase-space formulation. This can be seen by deriving the momen¬ 
tum variable 

dL 

which highlights the ambiguity between the position variables 6a and momentum variables TTa- 
The path integral over Grassmann variables therefore requires a slightly different approach to 
that over ordinary commuting variables. Different methods can be found, for example, in Berezin 
and Marinov [23], Swanson Junker |2B], and Aouachria |29| . 

Probably the simplest method of performing the Grassmann path integral is by using time- 
slicing. To achieve this, we will hrst express the Grassmann part of the Lagrangian (14.11) in 
terms of 6^, -ip = 6^ + id"^, and ip = 6^ — 

L{6^, 6^, Ip, ip, Ip, Ip) = iad^O^ + i-^'<P'4’ + + aB^{x, t)ipip, (4.16) 

where we have again assumed B(x, t) to be polarized in the z-direction, B(x, t) = (0,0, B^{x, t)). 
The Grassmann path integral that we now wish to evaluate is 

Gn jjj 'DO^V'ip'D'ip exp dtL(0^, 0^,-0, -i/i,-i/i)^ , (4.17) 

where Gn is an undetermined normalisation constant. 

The equations of motion for this Lagrangian are, of course, just those given in (14.8p and (|4.9I) 
and the integral over 6^ is trivial because of (|4.8p . Thus, the corresponding term effectively 
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vanishes. The effective Lagrangian now reduces to 


ip, ijj, Ip) = i^ipip + i^ipip + aB^{x, t)'ipip. 
Combining the two kinetic terms using an integration-by-parts, 

rT 


(4.18) 


S = 


= i 


f dtL{^p,^p,^p,^p) 

Jo 

a - - r -■ 

{'^T'lpT - Tpoipo) +ia dtipip - iB{x, 
^ Jo 


(4.19) 


where we have dropped the superscript on B{x,t), and where we can then simplify the action 
by a change of variables 


y = ^exp(z / dTil(x(T),r) ) =-0exp (i/l(t)). 


This action is now 


a 

-S' = iXTXT - XoXo) + / dtxx- 

^ Jo 


(4.20) 


(4.21) 


In the time-slicing method, the time-integral is replaced with a discretised version by splitting 
the time interval [0,T] into N slices of duration At = TjN such that 


pT N-i 

ia / dtxx ~ V Xn+i 
^0 n=0 


Xn+l Xn 

At 


At. 


(4.22) 


Expanding this summation, except for the end-points n = 0 and n = N, it is found that for a 
fixed index n = m there are three terms: 


( Xm-l-lXm + XmXm XmXm—l) 


(4.23) 


Using the usual properties of Grassmann variables and Berezin integration (Appendix [B]) , it is 
then straightforward to integrate the exponential of these terms over Xm and Xm] 


dXmdXm exp {ia {-Xm+lXm + XmXm “ XmXm-l)} = exp {-ZCCXm+lXm-l} ■ (4.24) 


Using this result for each time-point in the summation except the end-points, the N — 1 inte¬ 
grations over the pair Xn, Xn give 


{ia)^ ^exp{fa(xTXT - XTXo)}, 


(4.25) 


where the subscript N has been replaced by the continuous-time label T. Clearly, as A —>■ oo 
the term does not converge to a finite value and therefore some form of regularisation 

is required; this is the purpose of the constant Gn in (I4.17p . Setting Gat = GQ{ia)~^ and 
combining ()4.25l) with the boundary terms of the action in (|4.21l) . the path integral is 


— 2, 

p{xt,T I xo,0) = —Go exp ( -ia 
a 


XtXo - ipCTXT - ^XoXo 


(4.26) 
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Using (I4.20p to express p in terms of "0) we finally have 


pi^pT,T I V'OjO) 


—Go exp 
a 





-Tpri^T 



which agrees in form with Swanson’s evaluation m- 


(4.27) 


4.2 Fermionic density supermatrix 

Previously, when we calculated the purely bosonic path integral for the diffusion NMR experi¬ 
ment (section [231) , we interpreted the result G(xt, P; xq, to) as a classical propagator. However, 
we will now interpret the result of the fermionic path integral p{'il)T,T | ■ 00 , 0 ) as a statistical 
mechanical density matrix. 

Since p{'4)T-,T \ ?/^0)0) is constructed from the Grassmann elements ^|J and ip, the density 
matrix is actually a density supermatrix of even degree (see Appendix[C|). In the case of ordinary 
density matrices or functions of the form p {xt, T \ xo, 0), for xt, xq € M, the trace of p is taken 
by summing or integrating over all states xt such that xq = xt- Thus, 

p (T) = tr (p) = j dxTP (xt, T I XT, 0) . (4.28) 

For a supermatrix, the natural operation is not a trace but a supertrace (see Appendix lU]) . 
Since p is an even-degree supermatrix, the supertrace has the form of the difference of two traces; 

p = str (p) = tr(/9+)-tr(/3_), (4.29) 

where 

P+ = P (V’T, T I V'T, 0), (4.30) 

P- = p(V’r,T I -V’TjO) , (4.31) 


which is interpreted as the sum over all periodic, p+, and antiperiodic, p_, orbits. 
Calculating p± from (|4.27|) . 


P± = —Go exp 
a 






(4.32) 


which is similar to the state sum model of fermions on a circle with gauge group U(l) obtained 
by Barrett et al. [30]. Using the properties of Grassmann variables to expand the exponential. 


p± = — -Go — ioiipT^^T Tc ^ ^ 


(4.33) 


and therefore, 

P-i_ — p— = —2Go'0tV’t6 (4.34) 

The above expression for p_|_ — p_ is possibly related to the off-diagonal terms in the density 
matrix for a spin-^ particle. Since we have assumed that the magnetisation is purely in the 
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transverse plane, we could represent it by the Bloch vector u = (cos/3, sin/3,0). The density 
matrix in this two-state representation (| -|-), | —)) is 

1 1 / 1 

p=-(l + u.„) = -(^_,, 'J, (4.35) 

where a are the Pauli matrices. However, whether or not this is merely a superficial resemblance 
is unclear. 

Returning to the density matrix in (14.34^ . we compute the supertrace by integrating over 
the Grassmann variables; 

p{T) = str ip) = j j dV'rd^T {p+ - p-) = (4.36) 

from which it would be convenient to set Go = —1/2. 


4.3 Effective bosonic Lagrangian 


We began, in section 01 by proposing that the Lagrangian for the diffusion NMR experiment can 
be formed by combining the stochastic diffusion Lagrangian with the Grassmann Lagrangian, as 
in equation ()4.1h . and showed that the combined Lagrangian reproduces something close to the 
correct equations of motion. We then showed that it was possible to perform the path integral 
over the Grassmann variables, without placing any restrictions on the form of B{x,t) except 
that it was polarized in a single direction, which we took to be the ^;-direction. We are now in 
a position to recombine the result of the Grassmann path integration, equation (I4.36D . with the 
remaining bosonic path integral over The path integral is now 


Vx^Vx^Vx^ exp ( / dt — -^XaX°‘ ) p{T) 


1 


AD 



1 


Vx^Vx^Vx'^ exp [I dt — ) exp {—ij5{T )), 


(4.37) 


and using (|4.2UI) to express /3(T) in its integral form, it is easy to see that the effective bosonic 
Lagrangian is given by 

1 Q 

L{x,x) =--^XaX^ — iB^{x,t)- (4.38) 

This is precisely the bosonic Lagrangian given in equation (I2.15p of section 12.21 with the special 
form B^{x,t) = gait)x°-. 

Furthermore, we note that the actual form of the original, purely bosonic, Lagrangian— 
the real part of the Lagrangian—was not important. This means that more general bosonic 
Lagrangians, such as those with potentials, relaxation terms, and velocity fields are also possible. 
For example, 

L{x,x) = (x - v(x,t))a (x - v{x,t)T - U{x,t) - iB^{x,t), (4.39) 


where 1)“ and U are effective velocity and relaxation (or potential) terms which are defined 
in (|A.3[) and (IA.4p of Appendix O 
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5 Discussion 


Although we have mostly concentrated on the simplest case—diffusion in a linear field gradient— 
we have shown that the method of path integration over ordinary commuting position variables 
(as was demonstrated in section [2]) is appropriate to model the magnetisation of a diffusion NMR 
experiment. Except for some generalisations which were not considered, such as non-constant 
diffusion tensors, the only uncertainty in this approach is in the observation that the classical 
trajectories for position variables become complex-valued. This led to the formation of 

Conjecture 12. II which specified that the end-points of the trajectory should be real numbers. 

The presence of the complex trajectories suggested that there were some internal degrees of 
freedom which were associated with the processing magnetisation of an NMR experiment, and 
to model these degrees of freedom explicitly we used the anticommuting Grassmann variables. 
Grassmann variables have previously been used many times to model spin at the classical level, 
with some of the earliest work being by Berezin and Marinov [221124] . and Casalbuoni [23], and 
can be considered to be the classical analogue of spin as h ^ 0. There are other ways to model 
spin or processing magnetisation (such as with representation variables of SU(2) [25] . or as 
geometric algebra m. or based on the spherical top [32|) but here, despite the unfamiliarity of 
such objects in NMR, we have employed Grassmann variables because of their relative simplicity 
and their connection with other fields of physics and mathematics. 

By combining the commuting variables with the Grassmann variables 0*, and thus pro¬ 
ducing a pseudoclassical model, it was shown that parts of the imaginary component of the 
trajectory for x“(t) appeared to be of a purely internal nature. That is, the imaginary trajec¬ 
tory of the form A“(t) = Aq -|- has no effect on the trajectories of the Grassmann variables 
or their path integral and only serves to shift the end-points of the position variables such that 
+ iX^rp. Since such a shift will, in general, affect observables (such as the average 
magnetisation) and there appears to be no obvious physical interpretation of the imaginary shift, 
we conclude that the end-points x^rp are purely real-valued, in agreement with Goniecture 12.11 

There are a number of directions in which this work could be extended. As already men¬ 
tioned, no attempt has been made to model the effects of boundary conditions such as reflective 
or absorbing walls, permeable membranes, etc. Unfortunately, even in fairly simple cases the 
propagators are generally difficult to solve analytically. Although the perturbation methods of 
path integration possibly offer some interesting ways of approaching these cases, techniques such 
as those of Mitra et al. [33] and Ghadirian et al. [3l] might be more profitable. 

In this paper, we modeled the transverse magnetisation only and implicitly assumed that any 
radiofrequency pulses merely modified the effective time-dependence of the magnetic field. In 
many cases, this simplification is sufficient but, for example, in the presence of inhomogeneous 
radiofrequency fields this approach might not be adequate. The work of Aouachria [29l|35] has 
shown that radiofrequency fields can be straightforwardly included as potentials in the Grass¬ 
mann variables and it would be interesting to explore the interaction between such potentials 
and the stochastic part of the Lagrangian. 

The pseudoclassical model involves a Lagrangian which is a combination of commuting vari¬ 
ables and anticommuting variables 0*. Translating this into the parlance of particle physics, 
we have bosonic fields x“(t) and fermionic fields 9^{t) on the one-dimensional manifold of time, 
and the natural question arises of whether there could exist a supersymmetry [2011261136] between 
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these fields. There are some good reasons to think that there could be. 

The Torrey-Bloch equation is somewhat similar to the Pauli equation for a spin-^ quantum 
particle in a magnetic field, which is known to exhibit a supersymmetry [8]. Whilst the Torrey- 
Bloch equation is not simply an analytical continuation (Wick rotation) of the Pauli equation, it 
is known that the diffusion or Fokker-Planck equation can also possess a supersymmetry |26[l37j . 
In the absence of a magnetic field, it is easy to show that the Lagrangian —XaX°‘ + i9a0°‘ has a 
trivial supersymmetry—trivial in the sense that the dynamics of the bosonic and fermionic fields 
are completely independent and therefore uninteresting. Adding the magnetic field so that the 
two sets of fields couple, as with the Lagrangian (14.ip . however, destroys the supersymmetry. 
This is because supersymmetry requires that there is also some form of additional potential 
which is related to the magnetic field. In the case of the Pauli equation, the supersymmetry 
is a minimal Af = 1 supersymmetry and the related potential is the vector potential of the 
magnetic field [8] . In the case of Witten-type [26l[38ll39] M = 2 supersymmetry (which we have 
not explored here) the relation is less clear in the present context and requires clarification. 

It is fairly certain that supersymmetric Lagrangians can be constructed for both Af = 1 and 
Af = 2 supersymmetry but it is not yet clear whether these models correspond to physically 
realisable systems in the context of diffusion NMR. It is also not yet clear whether such su¬ 
persymmetries have any theoretical or experimental consequences. This is the focus of current 
research. 

Regardless of whether any form of supersymmetric model can be constructed the results 
presented here—for the ordinary (bosonic) path integral or for the pseudoclassical model—form 
the basis on which path-integral methods can be developed for diffusion NMR experiments. 
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A Appendix: Lagrangian of the Fokker-Planck equation 

Here is provided a derivation of the Lagrangian corresponding to the Eokker-Planck equation. In 
the process, we shall also discover the correct form of the Legendre transform and its associated 
equations of motion. The derivation follows closely that given by Wio [B]. 

We begin with a fairly general Eokker-Planck equation of the form 

|[ = {D{x,t)f) - V (v(x,t)/) - Uix,t)f, (A.l) 

where / = f{x,t), D{x,t) is a scalar diffusion coefficient, v{x,t) is a velocity field, U{x,t) is a 
sink or source term. In the context of diffusion NMR, U{x,t) may include transverse relaxation 
and the imaginary NMR (Bloch) term i'yB{x,t), as in equations (j2.8p or ()2.inh . Although this 
Eokker-Planck equation is quite general, for simplicity, we do not consider the more general 
cases of tensor diffusion or the effects of topological constraints (such as impermeable barriers). 
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Actually, the case of a constant diffusion tensor is straightforward but becomes more complicated 
if it is spatially or time-dependent. 

Expanding the derivatives in (lA.lh . we can rewrite it as 


where 


^ = D{x,t)S/‘^f - v(x,t) ■ V/ - U{x,t)f, 


V = V — 2VD, 

U = U + V - v -V^D. 


(A.2) 


(A.3) 

(A.4) 


Let f{x 2 ,t) be a solution of (IA.2D at x = X 2 . Then the solution at an inhnitesimal time 
increment e later is 

f{x 2 ,t + e) Ki f{x 2 , t) + e—J{x 2 , t). (A.5) 


With this, we can write, 


d 


f{x2,t + e) = ( 1 + e—)/(x2,t) 


j d^xi - X2)f{xi,t), 


(A.6) 


where a three-dimensional delta-function integral has been inserted to relate the position X 2 at 
time t -I- e to xi at time t. 

Using the Fokker-Planck equation (IA.2|) . we can replace the time-derivative in (IA.6h with a 
derivative operator which operates on X 2 . That is, 


/' 


f{x 2 ,t + €)= d^xi{l + e D{x 2 ,t)S/^ - v{x 2 ,t)-V - U{x 2 ,t) ) 6^{xi - X 2 )f{xi,t), (A.7) 


where V operates on X 2 only. 

Replacing the delta-function with its Fourier representation, 

- 2^2) = / (A- 8 ) 

we can effectively replace V with —ip in ()A.7p such that it becomes the phase-space integral 




D{x 2 ,t){-ipf - v(x 2 , t) • (-ip) - u{X2,t) 


„ip-{xi-X2 




Approximating the term in square brackets as an exponential, we can write this as 

(xi - X2) 




ip 


+ D{x 2 ,t){-ipY - v{x 2 ,t) ■ (-ip) - U {X2,t) 


(A.9) 

(A.IO) 
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We can now interpret the term in square brackets in (jA.lOp as the infinitesimal action integral 




S = 


rt+e 

dTL{x,x) = / dr 


—ip • X + t)(— — v(x, r) • (— ip) — C/(x, r) , (A. 11) 


where the term (x 2 — xi)/e has been identified as the velocity x. Prom this it is easy to see that 
the Lagrangian can be inferred to be 


L{x,x) = (—ip) • X — —D{x,t){—ip)'^ + v(x,t) • (—ip) + U{x,t) 
= {-ip) • X - i^(x, i-ip)), 


(A.12) 


where, since the first line has been deliberately written in the form of the Legendre transform, 
we can identify the Hamiltonian in the second line. Note that the Legendre transform above is 
equivalent to the usual form with p i— >■ —ip. From this, we can deduce that the Hamiltonian 
and Lagrangian equations of motion are 


dH 

dH 


= -IX 

dpa 

dL 

dL 


dx<^ ~ 


(A.13) 

(A.14) 


To complete the derivation so that L is expressed in 
Hamilton’s equations above to obtain 


terms of x and x only, we can use 


dH 

dpa 


-ix“ = 2D{x,t)p°‘ - ii}“. 


(A.15) 


and therefore, 


ip“ = 


V — X 

2D 


Substituting this into the Lagrangian ()A.12p we obtain the final result 


(A.16) 


L{x^ x) = 


(x — v(x, t))^ 
4D(x, t) 


U{x, t). 


(A.17) 


B Appendix: Berezin integration 

A brief introduction to the real Grassmann algebra was given in section 13.11 Here, we equally 
briefly explain integration over Grassmann variables, known as Berezin integration |20ll24lH0] . 

Let rj G r(l) be a real Grassmann variable. Then, since rf = 0, any function of rj, f{rj), can 
be expanded very simply as 

f{r])=a + hip (B.l) 

for constants a, 6, which may in general be of even or odd degree. To perform an integral over 
r], we need the differential dr/. This is also an odd Grassmann element and thus 

r/dr/= —dr/r/. (B-2) 
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This means that in defining the integral, a convention must be established. Here we use the 
following; 


j drj = 0 , 
dijT] = 1 . 


(B.3) 

(B.4) 


Note that Berezin integration is always indefinite. Applying the above to a function of t], we 
have 

J dr]f{rj)= J dr? (a + 5r?) = (— J dr/r? = ( —(B.5) 

Extending these rules to integration involving more than one Grassmann variable is straight¬ 
forward except that care must be taken with regard to the signs. Thus, for rji,ri 2 £ r(2), 


(B. 6 ) 
(B.7) 

(B. 8 ) 

(B.9) 
(B.IO) 

regardless of the degree of ci 2 - 

For more subtleties, in particular for complex Grassmann variables, see, for example, Swan¬ 
son [TT] . 


J dr]i = 0 , 

drjiTjj = 6ij, 

d?7i?7id?72r72 ^ ~ fj d?/idr72??i% 

dr7id772r72??i = 1- 

For a function of two Grassmann variables /(r/i,r? 2 ), which has the general expansion 

= a + + hm + ci2??i??2, 

d?7idr?2/(m)%) = / / dr/id?72Ci2??i% = -ci2, 


C Appendix: Supermatrices and supertraces 

Definition C.l (Supermatrix |21jb A matrix M is called a real (or complex) supermatrix if 
its entries have values in a real (or complex) Grassmann algebra T = Tg 0 Tj. More precisely, 
consider the set of [m + n) x [p + q) supermatrices M of the form 

"=(c d ) <C'1) 

where A, B, C, D are mxp, mxq,nxp and n x q submatrices respectively. The supermatrix 
M is called even (or of degree d) if A,D £Tq and B,C £ Tj, while it is called odd (or of degree 
1) if A, D £Ti and H, C € Tg. 


22 


Definition C.2 (Supertrace |21]i. The supertrace of the supermatrix M is defined by 


str(M) = tr(A) - (-l)'i®g^tr(D) 


tr(A)-tr(D) if Mg To 
tr(A)+tr(D) ifMGTi. 


(C.2) 
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